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Abstract 

We compute the thermal conductivity of color-flavor locked (CFL) quark matter. At temperatures 
below the scale set by the gap in the quark spectrum, transport properties are determined by collective 
modes. In this work we focus on the contribution from the lightest modes, the superfluid phonon 
and the massive neutral kaon. The calculation is done in the framework of kinetic theory, using 
variational solutions of the linearized Boltzmann equation. We find that the thermal conductivity 
due to phonons is ~ 1.04 x 10 26 ^ 00 Ag 6 erg cm -1 s" 1 K _1 and the contribution of kaons is 

k( k ) ~ 2.81 x 10 21 f£ 100 TjJ^y m 10 5//2 erg cm -1 s _1 K _1 . These values are smaller than previous estimates, 
but still much larger than (in the case of phonons) or similar to (for kaons) the corresponding values in 
nuclear matter. From the phonon thermal conductivity we estimate that a CFL quark matter core of a 
compact star becomes isothermal on a timescale of a few seconds. 



I. INTRODUCTION 



In this paper, we explore the thermal conductivity of quark matter in the regime of large baryon 
density and low temperature. It is expected that cold dense quark matter is a color superconductor, 
and that at asymptotically high density the ground state of three-flavor quark matter is the color- 
flavor locked (CFL) phase [1]. The work described in this paper is part of an ongoing research 
effort with the goal of determining the transport properties of not only the CFL phase, but also 
of other, less dense, phases of quark matter. Previous work has focused on the shear viscosity [2], 
bulk viscosity [3, 4], and neutrino emissivity of the CFL phase [5, 6]. A general discussion of the 
kinetics of a CFL superfluid is given in [7]. There are also calculations of the transport properties 
of unpaired quark matter [8-10], and some results for the transport properties of a kaon condensed 
CFL phase [11]. The long term goal is to connect these calculations of transport properties to 
possible observational signatures of high density quark matter phases in the core of compact stars. 

The thermal conductivity of dense matter plays an important part in the cooling of compact 
stars [12-14]. The cooling history of a star depends on the rate of energy loss by neutrino emission 
from the bulk and, at very late times, photon emission from the surface, on the specific heat, and 
on the thermal conductivity. Thermal conductivity determines how quickly different layers of the 
star become isothermal. It is generally assumed that the core of the star becomes isothermal very 
quickly, but that the outer crust is a poor conductor of heat and may take several hundreds of 
years to reach the temperature of the inner core [15-17]. In order to verify that this is true for CFL 
quark matter core we have to compute the thermal conductivity of CFL matter. There is an earlier 
calculation of the thermal conductivity of the CFL phase [18]. This calculation was based on a 
simple mean free path estimate. In the present paper we will perform a more definitive calculation 
based on the linearized Boltzmann equation. We will show that while the mean free path estimate 
is not reliable, the main conclusion of [18] is valid - a CFL quark matter core becomes isothermal 
on a very short time scale. 

The thermal conductivity, like other transport properties, depends on the properties of quasi- 
particle excitations. In the CFL phase up, down, and strange quarks are all gapped, with the gaps 
typically much larger than the appropriate temperature of the compact stars. Hence, the quarks are 
unlikely to have any significant contribution to transport properties. Below the gap the excitation 
spectrum consists of collective modes associated with the spontaneously broken symmetries of the 
CFL phase, see [19-22]. There is a phonon mode related to superfluidity, which we can view as 
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the Goldstone boson associated with the breaking of the U(l) of baryon number, and there is a 
meson octet coming from the breaking of the SU(3) C x SU(3)l x SU(3)r color and chiral flavor 
symmetry to the diagonal subgroup SU(3)c+f- This diagonal subgroup is the residual symmetry 
of the color-flavor locked diquark condensate. The quantum numbers of the meson octet coincide 
with those arising from the breaking of chiral symmetry in the QCD vacuum, and we will refer 
to these modes as pions, kaons, etc. However, the mass hierarchy of the meson octet in CFL is 
different, [20], such that the kaons are the lightest excitations. In this paper, we will calculate the 
thermal conductivity due to the lightest excitations in the CFL phase, phonons and kaons. 

This paper is organized as follows. In Section II and III we will outline the calculation of the 
thermal conductivity in kinetic theory, in Section IV we will discuss the interaction among the light 
excitations and derive the associated collision terms, and in Section V we will present numerical 
results. We end with some conclusions in Section VI. 



II. THERMAL CONDUCTIVITY AND TRANSPORT THEORY 

The thermal conductivity is defined in hydrodynamics as 

q = -nVT (2.1) 

where q is the heat flow and k is the thermal conductivity. In kinetic theory, the heat flux can be 
written in terms of the quasi-particle distribution function as 

d 3 p 



(27r) sVWp (2-2) 

where v p = dE p /dp is the particle velocity and Sf p = f p — f p is the deviation of the distribution 
function from the form in local thermal equilibrium, given by 

fp = e (p^-ri/T _ I ■ (2.3) 

Here, is the local fluid velocity and /i is the chemical potential (which is zero in the case of 
phonons and equal to the hypercharge chemical potential in the case of the kaons). Since we are 
only interested in contributions to the stress-energy tensor arising from thermal gradients we can 
write Sfp as 

f°(l + f°) 

6f P = - JA + Jp) g(p)p-VT, (2.4) 
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where g(p) is a dimensionless function of the magnitude of the momentum. Inserting this into 
Eq. (2.2), we can then read off the heat flux as 

VT f d?p 



3T 3 

and the thermal conductivity as 

K= W*I J^fi^ + ftBtovpEtP (2.6) 
Kinetic theory also gives some constrains on the form that 5f p can take. Conservation of particle 
number, energy and momentum imply [25] 



dT5f p = J dTE p Sf p = J dT P 5f p = , (2.7) 

where dT = j^s- Note that the number of phonons is not explicitly conserved. However, we will 
only consider 2 «-> 2 processes and the particle number constraint is expected to hold. We can show 
that for Sf p of the form given in Eq. (2.4) the first and second constraints (energy and momentum 
conservation) are automatically satisfied, but the third one (particle number conservation) is non- 
trivial. We get 

= J dT P 5f p = ^f 3 J dTf p (l + f p )g(p)p 2 . (2.8) 

Note that if v p is independent of p, then the contribution to the thermal conductivity vanishes 
due to the constraint. This implies that phonons with an exactly linear dispersion relation do 
not contribute to the thermal conductivity. This result was first derived in connection with the 
transport properties of superfluid helium [26]. The thermal conductivity of superfluid helium is 
dominated by rotons and phonon-roton scattering. Note that the phonon and roton are part of 
the same excitation curve. The name roton refers to the part of the excitation curve in which the 
disperion relation is very non-linear and E(p) develops a second minimum. In the following we will 
compute the thermal conductivity of the CFL phase due to non-linearities in the phonon dispersion 
relation, and due to massive kaons. 

To solve for the thermal conductivity, we need to determine the form of g(p). The quasi-particle 
distribution functions satisfies the Boltzmann equation, 

dfp _ 9f p df p df p 

ir-^t +Vp -^ + Fext ' - c[a (2 - 9) 

where F ext is an external force and C[f] is the collision term. In the case of binary scattering the 
collision integral is given by 



C[f P ] = ^r I W%P + K-K'- P')\M\ 2 D 2 ^ 2 , (2.10) 

ltj P Jk,k',p' 



where 

d 3 q 



, q , (2^2E q (2 - U) 
is the one-particle phase space integral, and D 2 ^ 2 = /p/jt(l + /p'X 1 + fk>) ~ (1 + + fk)f P 'fk' is 
the term involving the distribution functions for 2^2 scattering. Linearizing in 5f we can write 
Z?2^2 = -Do + SD where 



5D = _W± t*J± J JL± Ag .\j T ^ (2.12) 



T 3 
and 

A, = 0(p)p + </(£)*; - ^A;')*/ - <?(?>' • (2-13) 

Note that the collision integral will vanish for D = Dq. This is a consequence of detailed balance, 

which implies that \M.\ 2 is invariant under (p, k) <-> (p f k') and that fpf k (l + /£/)(! + fp') 
: 1 + /°)(1 + fk)fk'fp'- We can then write the collision integral as 

F • VT 

C[/ P ] = (2.14) 



where 



and 



fpfkO- + /$)(! + fp') A 
dr kik , jpl z— A g (2.15) 



J T 



[drk,k>j= f {2ir)%P + K-P' -K')\M\ 2 . (2.16) 

J Jk,k',p' 

We can simplify the left-hand side of the Boltzmann equation using our definition of the distribution 



function in Eq. (2.3). In general, we can write 

dj_ _ ■/?(! + /,° 
dt~ Qp T 2 



up P " , P V VT, (2.17) 



where the form of a p depends on the properties of the quasi-particles we are studying. In the 
following we will consider two cases. The first is phonons with a non-linear dispersion relation 



(2 2 \ 
1 + e^f-J, (2.18) 



where v = l/\/3 is the speed of sound, 7 controls the curvature of the dispersion relation, and 
e = 7T 2 /?; 2 is a dimensionless measure of the non-linearity. The parameter 7 was computed by 
Zarembo [27], who finds 7 = —11/(540 A 2 ), where A is the gap in the fermion spectrum. The 
second case is massive kaons with 
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K K) = \/ v kp 2 + m «^^ +mK ' ( 2 - 19 ) 



where vk = v. We will specify the kaon mass nix in Section IV B. The quasi-particle velocities, 

dE p 

— i cxi c g,i v en uy 

,2^2 \ „,2„ 



8E P ■ , 



r= „( 1 + 3£ ^j, „f>=^. (2.20) 



The explicit form of a p is derived in Appendix A. We find 



(^=4et,M^-— "l ^ = ^-^ (2 21) 

" V 7 y ' p 2m| 2m x ' 1 ] 



We have now expressed the LHS of the Boltzmann equation, Eq. (2.17), in terms of the function a p . 
The RHS of the Boltzmann equation, Eq. (2.14), is given in terms of the unknown function g(p). 
Once g{p) is determined the thermal conductivity is given by Eq. (2.6). In practice we will solve 
the Boltzmann equation using a variational procedure. This procedure is based on an expression 
for the thermal conductivity in terms of a suitable integral over the collision term, which we will 
now derive. Consider the following integral over the LHS of the Boltzmann equation, 

T I dT ft MP = J dT + f> 2a ^{p), (2.22) 

where c is a constant. We observe that the RHS of this equation has the same structure as Eq. (2.5) 
for the heat flux. For the two expressions to be equal, we need 

c J dT f p (l + f° p )p 2 a p g(p) = j dT f p (l + /> 2 ^ g(p) . (2.23) 

The two sides of this equation can be matched by using the constraint Eq. (2.8). The constraint 
implies that constant terms in v p E p /p and a v do not contribute to Eq. (2.23). We can then 
determine c by matching the p 2 terms in v p E p /p and a p . We find c = 1 for both phonons and 
kaons. Using this fact together with the Boltzmann equation we can express the heat current in 
terms of the collision integral 

J dTg(p)p-F. (2.24) 



VT 

Q = - 



3T 3 

The corresponding expression for the thermal conductivity is 



12T 4 

where 



k = 77tU / dr pAk ,, p ,Al. (2.25) 



/ dT pAk , p , = j (2tt) 4 5(P + K-K'- P')\M\ 2 f° p f° k (l + #)(1 + ft) (2-26) 

J J pkp'k' 

and we have used Eq. (2.15) and the symmetries of the matrix elements to derive Eq. (2.25). We 
now have two expressions for the thermal conductivity, Eq. (2.6) and Eq. (2.25). Both equations 



depend on the unknown function g(p), but the equivalence of the two equations rests on the fact 
that g(p) satisfies the Boltzmann equation. As we will now show, this fact can be used to determine 
gip)- 



III. VARIATIONAL SOLUTION TO THE BOLTZMANN EQUATION 

We will solve the Boltzmann equation by expanding g{p) in a basis of orthogonal polynomials. 
The procedure is variational in the sense that the result for n obtained in a truncated basis provides 
an upper bound on k for the exact solution. The expansion has the form 

g(p) = J2 b sB s (p 2 )i (3.1) 

s 

where B s (p 2 ) is a polynomial in p 2 of order s. The coefficient of the highest power is set to one. 
This means that B = 1, B\ = p 2 + c 10 , etc. The polynomials B s (p 2 ) are orthogonal with regard to 
the inner product 

J dVf p {\ + f p )p 2 B s {p 2 )B t {p 2 ) ee A s S st . (3.2) 

The functions B s {p 2 ) are a generalization of the Laguerre polynomials used in solutions of the 
linearized Boltzmann equation in classical physics [25]. Starting with B = 1, we can solve for all 
higher polynomials and their normalizations, A s . This is laid out in more detail in Appendix B. 
Inserting the expansion Eq. (3.1) into the constraint Eq. (2.8) we get 

= J dTf p {l + Op 2 B s (p 2 )B (p 2 ) , (3.3) 

where we have used Bq = 1. We conclude that the constraint is satisfied if bo = 0. Using the 
polynomial expansion in the first expression for k, Eq. (2.6), gives 

K = w^ hs f dTf " {1 + f0p) p2 Bs{p2) M + aiBl) = W* hlAl ' (3 - 4) 

where we have written v p E p /p = ao-Bo + ai-Bi- The coefficients ao and a\ are determined by 
Eq. (2.20). We get a\ = 4ew 4 /T 2 for phonons and a\ = v 4: /2m 2 K for kaons. Substituting the 
polynomial expansion into the second expression for the thermal conductivity, Eq. (2.25), gives 
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where 



M, 



St 



j dr p> k,y,p'Q s Qt, Q s = B s (p 2 )p + B s (k 2 )k - B s (k' 2 )k' - B s (p' 2 )p' . (3.6) 



Requiring Eq. (3.4) and Eq. (3.5) to be equal gives an equation for b s 



3T 3 



(3.7) 



which is equivalent to the linear equation 



Y,M st b t = ^-A l 5 sl . 



(3.8) 



This equation can be solved by inverting matrix M. We get 

M 



V 



A 1 M~ 1 







(3.9) 



This equation determines b\ and, using (3.4), the thermal conductivity. We find 



9T 2 



A\M^\ 



(3.10) 



where M^ 1 is the (1, l)-element of the matrix M, the constants a± are given below Eq. (3.4) and 
Ai is given in Appendix A. M is an infinite matrix, and in practice we solve for k by restricting 
the dimension of the matrix to a finite number N. This procedure can be shown to be variational 
in the sense that [28] 



K > 



4a? 



(Mi)' 



- 9T V £^AM st 

for any value of N. This right hand side can be written as 



(3.11; 



(3.12) 



where M is the the inverse of the truncated N x N matrix. The bound is saturated as N — > oo. 



IV. COLLISION TERMS 



The matrix M depends on the 2 <-> 2 scattering amplitudes. In this Section we will compute the 
phonon and kaon scattering amplitude using a low energy effective lagrangian for the CFL phase. 
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FIG. 1: Lowest order Feynman diagrams for the 2-phonon scattering amplitude. 



A. Superfluid Phonon 



The effective lagrangian for the phonon field <p is given by [29] 



7T 



: (<9 M W) 2 + 



(4-1) 



9^ — " T T/ 108/^ 

where v — 1/ a/3 is the speed of sound and fi q is the quark chemical potential. We have displayed 
the leading three and four-phonon vertices. Higher order terms include higher powers of <ft or 
additional derivatives. These terms are suppressed by powers of the typical momentum over the 
quark chemical potential. The speed of sound as well as the coefficients of the three and four-phonon 
vertices are given to leading order in the strong coupling constant. The lowest order diagrams that 
contribute to phono n-phonon scattering are shown in Figure 1. 

The corresponding matrix elements were previously computed in connection with the phonon 
contribution to the shear viscosity of the CFL phase [2]. The matrix element A4(P, K\ P'K') where 
P, K and P', K' are the in and out-going four momenta can be written as the sum of the contact 
term Ai c and the s, t, u-channel phonon exchange diagrams ftA s ,t,u- The individual terms are 

iM c = A [(P ■ K){P' -K) + (P- K'){P' -K) + (P- P'){K ■ K')} 

iM s = g 2 [2(p + k )P ■ K + Po K 2 + k P 2 ] [2(p' + k' )P' ■ K' + p' K' 2 + k' P' 2 ] G(P + K) 
iM t = iM s (P^K') 

iM u = iM t (P <-> K), (4.2) 

where G(Q) is the phonon propagator and the coupling constants A and g are given by 

4tt 2 2tt 



A 



108/4 
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(4.3) 



~q "t~q 

At leading order the phonon propagator is given by G(Q) = (q 2 — E 2 )~ l with E g = vq. The thermal 
conductivity is sensitive to non-linearities in the dispersion relation and we will use E q as given in 



S 



Eq. (2.18). We note that non-linear effects are generated by higher derivative terms in the effective 
lagrangian. Except for the non-linearities in the dispersion relation the structure of these terms 
has not been determined. We have verified that the matrix elements M st are not sensitive to higher 
order corrections in the vertices. We also note that a non-linear dispersion (with e < 0) regulates a 
possible on-shell divergence of the matrix element, and as consequence there is no need to include 
self-energy corrections in the fermion propagator as in [2, 30]. 

We can now insert the scattering amplitude into Eq. (3.6) and compute the matrix elements 
M st . Using the energy and momentum conserving delta functions as well as spherical symmetry 
the 12- dimensional phase space integral in Eq. (3.6) can be reduced to a 5-dimensional integral. 
First, we use the momentum conserving delta function to integrate over p'. We can choose p to lie 
along the ^-direction and perform the corresponding angular integrals. The energy conserving part 
of the delta function can be used to fix the magnitude of k' . Finally, we note that the integrand only 
depends on the relative azimuthal angle of k and k' and we can trivially integrate over the other 
azimuthal angle. The remaining 5-dimensional integral is computed numerically using standard 
methods, such as the VEGAS routine [31]. 



B. Massive Kaon 



The effective lagrangian for the meson octet is given by [19-21] 

£ = i^Ti [VoSVoS f -w 2 V£V£ f ] + ^Tr [M -1 |M|(E + E f )] (4.4) 

where £ = exp(i9 a t a / f n ) is the chiral field, t a with a = 1,...,8 are the SU(3) Gell-Mann 
matrices, f n is the pion decay constant, M = diag(m u , ma, rn s ) is the quark mass matrix, and 
\M\ is the determinant of the mass matrix. The covariant derivative of the chiral field is given by 
VoS = BqII— i[A, E] where A is the effective chemical potential A = —M 2 /(2n q ). At asymptotically 
high density the constants f n , v, and a can be determined by matching the effective theory to 
perturbative QCD [20, 21, 32] 

£2 21-81n2/i 2 1 3A 2 

where A is the fermionic energy gap at zero temperature. For m s ^> m u , the lightest excitations 
in the theory are kaons. The leading terms in the effective theory involving kaons are Cq + £j„t 
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FIG. 2: Spectrum of energies for the charged and neutral kaons. The mass splitting between K° and K + 
is given by 5m and the larger splitting between the kaons containing s and those containing s is given by 
m 2 s /2n q . 



with 



Co = l -{d^K°){d»K ) + l -{d tl K + ){d»K~ 



-m 2 K oK°K - -ml + K + K 



c 



lilt 



A (K°K°) 2 + A + (K + K~Y + (A + A+) (K°K°) (K + K~) 



where A . 



21 

m 2 K ^ + /fl and 



am u {m d + m s ) , m 2 K + = am d {m u + m s ) 



(4.6) 



(4.7) 



We have not included kaon interactions that involve derivatives of the kaon field. These terms arise 
from expanding the first term in (4.4), see [33]. For energies and momenta that are small compared 
to the kaon mass these terms are suppressed as compared to the terms in Eq. (4.6). 

The spectrum of kaons is shown schematically in Figure 2. The mass difference between the 
isospin doublets (K°, K + ) and (K°, K~) is proportional to m 2 / (2/x 9 ), while the mass splitting with 
the doublets is 5m 2 = am s (md — m u ). In the following we will only include the lightest kaon, the 
K°. The relevant interaction term is 



m 



K° 



hit 



24/* 



and the corresponding scattering amplitude is 

\M\ 2 



m 



K° 



f 4 



(4.9) 
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This expression can be inserted into Eq. (3.6) in order to compute the matrix M st . As in the case 
of phonon scattering, the phase space integral can be reduced to a 5-dimensional integral over the 
magnitudes of p and k, the polar angles cos(#fc) and cos($k') describing the angles between p and 
the respective momenta, and one azimuthal angle. 



V. RESULTS 



Before we show numerical results for the thermal conductivity, we would like to discuss the 
dependence of n on the dimensionful parameters that appear in the problem, the temperature T, 
the quark chemical potential n q , the fermion gap A, and the mass of the kaon. We begin by studying 
the scaling of the matrix elements M st . In the phonon case, we can rescale all momenta by the 
temperature and write M st as a product of the dimensionful parameters and a dimensionless phase 
space integral. A subtlety arise due to the possible role of on-shell divergences. There is no on-shell 
sensitivity in the phonon contribution to the shear viscosity and the corresponding matrix elements 
exhibit the naive scaling behavior, but the phonon mean free path is sensitive to the phonon self 
energy near the on-shell point [2, 30]. We find that the thermal conductivity also depends on the 
behavior of the phonon propagator near the on-shell point. In the present calculation this behavior 
is governed by non-linearities in the dispersion relation. Consider the contribution to M st from 
nearly on-shell phonons in the s, t or w-channel. We find 

(P) rp2s+2t+U j< dq 2 dT' rp2s+2t+U j rp2s+2t+12^2 



J (q 2 + ef(q,T')y // |e| 1 ' ) 

where q represent the momentum transfer in the relevant channel, dT' is a phase space integral over 
the remaining momenta, and f(q, V) is a function of these variables. 

In the kaon case we can assume that the temperature is smaller than all the other parameters 



and write the dispersion as E q = + q / (2m k). We rescale the momenta by y/m^T and extract 
all dimensionful parameters from the phase space integral. We get 

M™ ~ ^L( mK T) s+t ^/ 2 e- 2&m l T ~ ^( mK T) s+t+9 ^e~ 25m / T (5.2) 

Using Eq. (5.1) and Eq. (5.2), as well as the scaling of the normalization factor A\ derived in 
Appendix B, we find 

2§ phonons, 

-*f- \ / — kaons. 
V mK 
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FIG. 3: Numerical results for the thermal conductivity of the CFL phase due to phonons and kaons as 
a function of the basis size N used in the variational procedure. Results are shown for two different 
temperatures T = 0.2 MeV and T = 1 MeV. We have used fj, = 400 MeV, A = 100 MeV for the phonons 
and rriK = 10 MeV, (j,k = 5 MeV and /„- = 100 MeV for the kaons. Numerical results are shown as stars 
and crosses, the lines are only serve to guide the eye. 



We note that the contribution to the thermal conductivity from phonons is independent of 
temperature. The naive power counting gives k ~ fi 8 /T 6 , but the contribution from nearly on- 
shell phonons modifies the naive scaling behavior by a factor e 3 . Two powers of e come from the 
a\ term and one power of e comes from the nonlinear dispersion cutting off the collinear regimes 
in the propagators. Using e ~ T 2 then leads to a thermal conductivity that is independent of 
temperature. 

The thermal conductivity due to kaons is not exponentially suppressed by exp(— m^/T) as 
one might naively expect, because the exponential factors in A\ and M n cancel at leading order. 
This result is analogous to the well known fact that the transport properties of a dilute gas are 
independent of the density of the gas. We should note, however, that in practice the range of 
validity of this result is limited to temperatures that are not too small. At very small temperature 
the kaon mean free path is bigger than the system size, and heat transport is no longer a diffusive 
process. 

Numerical results for the thermal conductivity are shown in Figs. 3 and 4. In Fig. 3 we show 
the dependence of k on the number of basis functions. We observe that the convergence in the case 
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T (MeV) 



FIG. 4: Numerical results for the contribution of phonons and kaons to the thermal conductivity of the 
CFL as a function of temperature. The results shown in this figure were obtained using the variational 
method with N = 5. We have used (j, = 400 MeV, A = 100 MeV for the phonons and ttik = 10 MeV, 
UK = 5 MeV and f w = 100 MeV for the kaons. 

of the kaon contribution is excellent - using N = 1,2 gives results within a few percent of the final 
result. Also, the rate of convergence is independent of the temperature. In the case of the phonon 
contribution the rate of convergence is slower, but going up to N = 5 appears to give converged 
results. 

Fig. 4 shows the thermal conductivity as a function of temperature. We have used the results 
for N = 5. We observe that thermal conductivity due to kaons follows the square root dependence 
predicted by Eq. (5.3). Also, the thermal conductivity due to phonons is approximately temperature 
independent, in agreement with Eq. (5.3). The numerical results are well represented by 



A" m K V % 

We observe that for trk in the few MeV range the phonon contribution is significantly larger than 
the kaon contribution. 





(5.4) 
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VI. CONCLUSIONS 



We can convert these result into CGS units. We get 

> 1.04 x 10 26 tj™ , * 2.81 x 10 21 /^oo , (6-1) 

Ag cmsK ' cmsK 

where is the quark chemical potential in units of 500 MeV, A 50 is the gap in units of 50 MeV, 
f n ,wo is the pion decay constant in units of 100 MeV, Tm c v is the temperatures in units of MeV 
and 77iio is the kaon mass in units of 10 MeV. 

These numbers can be compared to a number of results in the literature. The thermal 
conductivity of the CFL phase was previously studied by Shovkovy and Ellis [18]. Based on 
a mean free path estimate they find that the contribution due to phonons is k,^ ~ 0.7 x 
10 32 T^ eV i?o,km cl ^ r & g K , where -Ro,km is the phonon mean free path in units of 1 km. This particular 
number was chosen because the authors argued that the mean free path is so long that it is effectively 
cut off by the size of the quark matter core. In the present work we have demonstrated that in 
the framework of kinetic theory there is no simple connection between the phonon mean free path 
and the thermal conductivity. We note, however, that if one replaces the mean free path in the 
Shovkovy- Ellis estimate with a more accurate result from kinetic theory, which is on the order of 
millimeters rather than kilometers [2], one obtains a thermal conductivity which is close to our 
result. Shovkovy and Ellis also argued that there is a large contribution to k from photons. The 
photon mean free path is indeed on the order of the size of the star [5]. This implies, however, that 
photons are not thermally coupled to the a CFL quark matter core. 

The thermal conductivity of unpaired quark matter at low temperature is k = 0.5m|,/a 2 [9], 
where mo is the gluon screening mass and a s is the strong coupling constant. Using mo — 500 
MeV and a s = 0.3 this gives k ~ 10 4 MeV 2 , about 5 orders of magnitude larger than the thermal 
conductivity of the CFL phase. The thermal conductivity of nuclear matter was studied in [34-38] . 
A typical value at nuclear matter saturation density and T ~ 1 MeV is k ~ (10 20 — 10 21 ) ei ^ r & g K , 
which is about 5 orders of magnitude smaller than the thermal conductivity of the CFL phase. 

We can convert the thermal conductivity to a timescale for thermal diffusion. The thermal 
diffusion constant is \ = K / C v, where cy is the specific heat. The time scale for thermal diffusion 
over a distance R is 

c vR 2 , a N 

r ~ . (6.2) 

K 

Using c v = 2vr 2 T 3 /(15w 3 ) for phonons, as well as T = 1 MeV, fi = 500 MeV, A = 50 MeV, and 
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R = 1 km we get r ~ 10 s, indicating that a CFL quark core will become isothermal rapidly. 

In this paper we have studied the thermal conductivity of the CFL phase using a linearized 
Boltzmann equation including a collision term that involves 2 <-> 2 scattering between phonons or 
kaons only. We have not studied a coupled transport equation for phonon-kaon scattering. Given 
that the relaxation times of phonons and kaons are quite different, it should be possible to find 
approximate solutions to the coupled problem. We have only provided a simple estimate for the 
relaxation time of temperature gradients. More accurate estimates will require a detailed model 
of the initial temperature profile, and a full solution of the dissipative two-fluid hydrodynamic 
equations for a CFL superfluid. 

Acknowledgments: We would like to thank C. Manuel and M. Manarelli for very helpful 
communications. This work was supported in parts by the US Department of Energy grant DE- 
FG02-03ER41260. 



APPENDIX A: STREAMING TERMS IN THE BOLTZMANN EQUATION 
1. Superfluid Phonon 

In this appendix we compute the quantity a p , which is related to the streaming terms (the 
left-hand side) of the Boltzmann equation. We will use the dispersion relation for the phonon given 
by 

(2 2 \ 
1 + e^j . (A.l) 

The left-hand side of the Boltzmann equation can be written as 

f = f + lV V/ p , (A,) 

where v p is the particle velocity. In local thermal equilibrium the time dependence of f p arises from 
the time dependence of the local temperature and fluid velocity. We have 

df P _ J° p (l + f° p ) (E ]L dT +p du\ ^ (A3) 



dt T \T dt dt J 

where we have used the fact that by going to the local rest frame we can set the local fluid velocity 
(but not its derivatives to zero). We have also used that the number of phonons is not conserved, 
and the phonon chemical potential is zero. The spatial derivatives can be written as 

v P ■ Vf p = - W + ffl (Et Vp .VT + v p . V(p ■ u)) . (A.4) 
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We can simplify these expressions using the Euler equation 

du t „ ldP„_ /A . 

— = — VP = — VT A.5 

dt p pdT K J 

where P is the pressure and p is the mass density. The mass density is defined by 

7T = pu , (A. 6) 

where 7r is the momentum density of the fluid. We can now collect all terms that contain a gradients 
of T. We find 

Noting that u p = {v p /p)p and using our definition of a p in Eq. (2.17), we can write 

T dP _ w„ . . . 

a p = -— -P p ^. A.8 
p aT p 

The pressure of the phonon gas is given by 

P = -T J dT In (1 - e -^ /T ) , (A.9) 

leading to 

rfP 1 wo 2tt 2 T 3 / 60e7r 2 \ 

= 3T 7 *W + ^ * IfoT i 1 + ~ J ' (A ' 10) 

where we have kept terms up to linear order in e. The mass density is given by 

p = ap- y ^ 2/ ° (i + # * w (i + 20e7r2) • (A - n) 

Inserting these results into Eq. (A.8) we get 



a< p > = 4et, 2 



vV 20vr 2 



(A.12) 



T 2 7 

We note that a p = for e = 0. This result is consistent with Eq. (2.8). In the case of phonons with 
an exactly linear dispersion relation, the LHS side (the streaming term) of the Boltzmann equation 
is zero. The RHS (the collision term) is zero if Sf p is a zero mode of the linearized collision operator. 
The zero mode associated with momentum conservation satisfies the constraints in Eq. (2.7) and, 
because of Eq. (2.8), gives a vanishing thermal conductivity. For e ^ both the LHS and the RHS 
are not zero, and the corresponding solution of the linearized Boltzmann equation leads to a finite 
thermal conductivity. 



1(3 



2. Massive Kaon 



In this section we compute the streaming term for kaons. There are two important differences as 
compared to the phonon case: kaons are massive, and they carry conserved charges (hypercharge 
and isospin) which couple to chemical potentials. For simplicity we consider only the neutral kaon, 
which couples to a single chemical potential (a linear combination of isopsin and hypercharge). We 
will also consider the limit of low temperature, such that T <C m, /a, 5m = m — u, where m = 
and /i = fix- The left-hand side of the Boltzmann equation can be written as 

dfv df p du df„dT df p da , . 

~dt = ^ • ~dt + df ~dt + ^Tt +V v gradient termS (A ' 13) 

The time derivatives of the temperature and the chemical potential can be converted to spatial 

derivatives using the the continuity equation. These terms do not contribute to the thermal 

conductivity. The derivative of the velocity field can be rewritten using the Euler equation 

du 1 

— = — VP. A.14 

dt p K 1 

Spatial derivatives of the velocity field only contribute to the shear viscosity. The spatial derivative 
of the chemical potential can be rewritten using the Gibbs-Duhem relation dP = ndfi + sdT, where 
n is the density and s is the entropy density. We get 

V/i = -VP--VT. (A.15) 

n n 

Collecting gradients of the temperature and pressure we find 

df P _ ( df p sdf p \ fv p df p ldf } 



dt \ dT n dfi J p ^ ^ \n dp: p du J ^ ( ^) 

At this point we need to specify explicit expressions for the thermodynamic quantities. The pressure 
of an ideal gas of massive bosons is given by 

P = -T J dT In (1 - e -(^V T ) . (A. 17) 

The entropy density, particle density, and mass density are 

s = /<*r/£(3(£p-/') + 'w)= + 

I Til 

P = ^ J dT P "f p (l + f p )^-n, (A.18) 

17 



where we have given analytical results in the low temperature limit. We also need the derivatives 
of the distribution function with respect to the thermodynamic quantities, 

dT T T 

df . . + f° P ) 



du T 

The coefficient of the VP term in Eq. (A. 16) is given by 



dp T 
Qt fOn + fO\ 



V 



P df ldf + {v t _m\ (A2Q) 



n dp p du T \n p J 

Using v p = v 2 p/E p ~ v 2 p/m we can see that v p /p ~ v 2 /m. We can also show that n/p = v 2 /m. 
This implies that the coefficient of VP is zero. The VT term is 

df {df sdf\_ „^ /°(1 + f°) (E p - n s 



-j- = — —-7T- ) Vp ■ VP = Jpy m Jp/ ^—^ — — ) v p - VP . (A.2i; 
dt \dT ndpj p T \ T n> r 



We can then read off the coefficient a p . We find 



a 



f_ \ J7I Ts] _ v 4 p 2 5v 2 T 



p m 



Ts 

E p - p 

n 



2m 2 2m 

APPENDIX B: ORTHOGONAL POLYNOMIALS 



(A.22) 



In this appendix we collect some explicit expressions for the orthogonal polynomials B s (p 2 ) 
introduced in section III. The starting point is Eq. (3.2) 

J dT7j(l + f p )p 2 B s {p 2 )B t {p 2 ) = A s S st , (B.l) 

together with Po = 1- We write 

Po = 1 

Pi = p 2 + C10 

P 2 = P 4 + c 22 p 2 + c 20 , (B.2) 

P 3 = p 6 + c 34 p 4 + c 32 p 2 + c 30 

and iteratively determine the coefficients c st . For example, we find 



20tt 2 T 2 
7T 



phonons, 



C10 = { /I r ~ (B.3) 
-^1 kaons. 
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which implies that B\ oc a p . This relation plays a role in ensuring the consistency of the linearized 
Boltzmann equation. The variational function g(p) = Bq is a zero mode of the linearized collision 
operator. We already showed that a term g(p) oc B is eliminated by the constraints, Eq. (2.8), 
and that it does not contribute to the thermal conductivity. The collision operator acting on 
any function g(p) is orthogonal to B , because B is a zero mode, and the collision operator is 
hermitean. This implies that the streaming term must be orthogonal to Bq also. The streaming 
term is proportional to a p , and orthogonality to Bq follows from the relation a p oc B\. 

The final result for the thermal conductivity requires the normalization constant A±. Using the 
result for B\ we get 

rT 9 phonons, 
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